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ABSTRACT 

It has been argued that blobs ejected from advection-dominated accretion flow through the accretion-ejection 
instability undergo expansion due to their high internal energy density. The expanding blobs interact with their 
surroundings and form strong shock, which accelerates a group of electrons to be relativistic. Then flares are 
formed. This model has advances in two aspects: shock acceleration and self-consistent injection. We derive an 
analytical formula of the injection function of relativistic electrons based on the Sedov's law. We calculate the 
time-dependent spectrum of relativistic electrons in such an expanding blob. The light-travel effect, the evolution 
of the electron spectrum due to energy loss, and the escape of relativistic electrons from the radiating region are 
considered, as well as the expansion (at sub-relativistic speed) of the coasting blob. A large number of light curves 
spanning wide spaces of parameters have been given in this paper. Regarding the symmetry, relative amplitude, 
duration of a flare, and the time lag between peak fluxes, we find four basic kinds of light curves for the non- 
expanding blob, and seven basic kinds of light curves for the expanding blob. The expansion weakens the magnetic 
field as well as enlarges the size of the blob. The predicted light curves thus show very complicated properties 
which are composed of the basic light curves, because the physical conditions are changing in the expanding blob. 
For the rapid decays of magnetic filed (for example, B oc r~"), we find the falling profile is a power law of time as 
uF v oc f( 1-a W 5 , where a is the power-law index of injected electrons. This decline is controlled by the decay of 
magnetic field rather than the energy losses of relativistic electrons. We also calculate the evolution of the photon 
spectrum from both non-expanding and expanding blobs. Different shapes in the phase of decreasing luminosity 
are then obtained for different parameter values. The photon index, a p h, keeps constant for non-expanding blobs 
when luminosity decreases, whereas a p h continues to decrease after the luminosity reaches its maximum for 
expanding blobs. It is expected that we can extract the information of ejected blobs from the observed light curves 
based on the present model. 

Subject headings: radiation mechanism-synchrotron-quasar-blazar 



Scott, & Vestrand 1978). The free expanding case was consid- 
ered under the assumptions that the electron spectrum is steady 
without time evolution (Jones & Tobin 1977; Vitello & Pacini 
1977 and 1978; Konigl 1978; Marscher 1980; Band & Grind- 
lay 1986). Then the effect of electron spectrum evolution with 
presumed injection function was generally considered (Celotti, 
Maraschi, & Treves 1991; Wang & Zhou 1995; Mastichiadis 
& Kirk 1997; Kirk, Rieger, & Mastichiadis 1998; Dermer, 
Sturner, & Schlickeiser 1997; Dermer 1998; Georganopoulos & 
Marscher 1998; Chiaberge & Ghisellini 1999; Li & Kusunose 
2000; Kusunose, Takahara, & Li 2000). The random fluctua- 
tions of the continuum from blazars are generally attributed to 
the propagation of relativistic shock (Marscher & Gear 1985). 
Large flares most probably result from a blob (Dermer, Sturner, 
& Schlickeiser 1997; Chiaberge & Ghisellini 1999) as consid- 
ered by the authors listed above. It should be cautioned that 
all the calculations avoid three important processes: 1) how to 
accelerate the electrons in the blob, 2) how to inject the accel- 
erated electrons into the blob if radiation regions are separated 
from acceleration regions, and 3) why and how to eject a blob. 
It is generally believed that these problems are related with ac- 
cretion disks around supper massive black holes as a central 
engine. The status of accretion disks may determine the prop- 
erties of the ejected blobs. 
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1. INTRODUCTION 

The time variations of received fluxes and emission spectra 
of blazar-type active galactic nuclei (AGNs) have received great 
attention from both theory and observation during the last 30 
years, because the multiwavelength observations of flares from 
blazars provide invaluable information how the central engines 
work (Bregman 1990; Urry & Padovani 1995; Utrich et al. 
1997). The flare variabilities from blazars are very common. 
However, the properties of these flares are diverse, and the out- 
bursts of the flares are over multiwavelength. Although we have 
a plenty of observational data of these phenomena, our under- 
standing of the mechanisms of the flares are still insufficient. 

Many attempts have been made to fit the observed continuum 
spectra by varieties of theoretical models. These calculations 
can be mainly divided into two kinds: homogeneous models 
(blob models) and inhomogeneous models (elongated jet mod- 
els). At early stage the homogeneous models were frequently 
invoked owing to their concise physical scenario. Rees (1966, 
1967) and Ozernoy & Sazonov (1968) originally studied the 
emergent spectrum from an expanding blob disregarding the 
evolution of the electron spectrum and blob dynamics. Sub- 
sequently, the expanding process was studied (Mathews 1971; 
Vitello & Salvati 1976; Canuto & Tsiang 1977; Christiansen, 
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The properties of the iron Ka line emission in radio-loud 
quasars may provide a new clue to understand the flare pro- 
cess, because it is a powerful probe to test the ionization state 
of matter close to black holes. The available data show that the 
line profiles in radio-loud quasars are very narrow (Sambruna 
et al. 1999; Reeves & Turner 2000), which are quite different 
from radio-quiet quasars and Seyfert galaxies. Thus an opti- 
cally thin advection-dominated accretion flow (ADAF, Narayan 
& Yi 1994) or ion-pressure supported tori (Rees et al. 1982) 
may power the central engine in radio-loud quasars, which may 
be explained by the properties of iron Ka line (Sambruna et 
al. 1999). If the accretion powers the central engine, then the 
accretion rate is roughly of order 



m ~ 10" 3 'V^L A7 M S , 



(1) 



where m is the dimensionless accretion rate normalized by the 
Eddington rate, V\q is the Doppler factor normalized by 10, 
i.e., 2?io = I?/10, L47 is the bolometric luminosity normalized 
by 10 47 erg/s, i.e., L 47 = L bol /10 47 , and M 8 is the black hole 
mass normalized by 10 8 Mq, i.e., M 8 = M B h/(10 8 Mq). Such 
a crude estimation is in agreement with that m <~ 10" 2 ~" 3 in 
FRII galaxies evaluated by Ghisellini & Celotti (2001). Indeed, 
such a low accretion rate implies that optically thin ADAF or 
ion-pressure supported tori may power the central engine. The 
observed bolometric luminosity suggests that the emission re- 
gions are highly relativistically boosted, which dominates over 
the emission from accretion disks (or ion-pressure supported 
tori). 

Tagger & Pellat (1999) showed there is a strong accretion- 
ejection instability, which increases with the height of an ac- 
cretion disk. Based on this instability, Wang et al. (2000) sug- 
gested that there is a large flare if a blob ejection takes place 
in optically thin ADAF of blazars. Comparing with the blob's 
surroundings, Wang et al. (2000) also suggested the observa- 
tional consequences of an ejected blob: nonthermal flares and 
recombination line emission, which are caused by the expan- 
sion of a blob owing to its high internal energy density. In fact 
there is an ample of evidence for the blob ejection from the nu- 
cleus. Especially, recent observations of 3C 273 by Chandra 
discovered the blob ejection from the central region (Marshall 
et al. 2001). If the flare arises from blob ejection, according to 
the model of Wang et al. (2000), what are the properties of the 
ejected blob in multiwavelength continuum? It is the goal of 
our paper to calculate the multiwavelength light curves of the 
emission from the ejected blob. 

This paper is arranged as follows. Section 2 describes our 
model and the formulation, and an analytical injection function 
is derived from Sedov's self-similar expansion. Extensive qual- 
itative analyses of light curves are given in §3. To compare with 
observations we provide a large number of light curves and con- 
firm the qualitative analysis in §4. Discussions and conclusions 
are given in the last section. 

2. THE DESCRIPTIONS OF MODELS 

Supposing that the inter-cloud medium is homogeneous, the 
expansion of a blob is quite similar to the 7-ray burst fire- 
ball. As shown by Wang et al. (2000) the expansion is sub- 
relativistic. The expansion leads to the formation of strong 
shock. Thus such an expansion can be described by the self- 
similar solution after the timescale t\ (Sedov 1969), which is 
rather shorter than that of the expansion phase. Sedov (1969) 



gives the self-similar solution: the shock speed v s h oc t~ 3 ^ 5 and 
shell radius r s h oc t 2 l 5 . Some electrons in the swept medium by 
blob expansion will be accelerated to be relativistic. Therefore, 
the expansion acts as the injection of relativistic electrons. We 
follow the usual assumption that a power-law distribution (7"") 
of electrons is formed via shock acceleration, where 7 is the 
Lorentz factor of an electron. The acceleration timescale is as- 
sumed to be much shorter than the expansion and synchrotron 
timescales (for the acceleration timescale, e.g., Kirk et al. 1994 
and references therein). Then the injection function Q(j, t) can 
be written as <2(7,f) ex r 2 h v s h <x t '/ 5 . Because the expansion of 
the blob will be drastically decelerated after long enough time, 
we introduce a cutoff factor e~'/' 2 , 



Qh,t) = Qor a t 1/5 e-'^., 



(2) 



where Qq can be obtained by the total injected energy Eq after 
integrating the time and electron energy, 



e ; 



6(2- a)E Q 



5m e c 2 [(7°ax) 2 - a -(7D 2 - Q K /5 ' 



for ail. (3) 



For a = 2, we have Q w 6£ , /[5ffJ e c 2 f 2 / ln(7 ma x/7min)]. Equa- 
tion (2) implies that the injection is roughly constant within a 
timescale t^- The evolution of electron energy can be deter- 
mined by 7 = 70/(1 +ajot), where 70 is the initial energy of 
electron and a = cajB 2 /6irm e c 2 = 1.29 x 10 _9 B 2 . Here c is the 
light speed, cr T is the Thomson cross section, m e is the electron 
mass, and B is the magnetic field strength in units of Gauss. 
The electron continuity equation in the homogeneous expand- 
ing blob reads 



at dj T esc 



(4) 



where b(j) = -aj 2 and T esc is the timescale of electron escape 
from the radiation region, which determines the later evolution 
of a flare showing up in lower frequencies. If re-acceleration 
of relativistic electrons takes place the pile-up distribution will 
be formed in the shell (Schlickeiser 1984), but we neglect the 
re-acceleration here. Therefore the solution of the continuity 
equation of electrons is given by (Ginzburg & Syrovatskii 1965) 
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7max 



^7,0= 77- / exp G(7o,f-r)^7o, (5) 

where r = b~ l (^')d^i' = (7" 1 -70 1 ) /a. The number density 

of the electrons reads n(j,t) = N(j, t) / (47rr 2 h Ar s h). These equa- 
tions describe the evolution of the relativistic electrons. 

It is known that the observed light curve is "contaminated" 
by two factors: the size of a blob and the evolution of the rela- 
tivistic electrons. Therefore the information of geometric struc- 
ture and motion is hidden in the light curves as well as the pure 
evolution of electrons. Figure 1 shows the geometry of the ex- 
panding blob in its coasting frame. Suppose at a given time ?biob 
(in the blob frame), the location of its center is C", and the dis- 
tance to the observer is d ' c ' . For the arbitrary point P(r sh , 9' , <j>') 
on the shell, the distance of P to the observer reads 



1- 



d , , 

o'c' 



COS?/> 



(6) 



where cos?/> = sin#'sin(0' + 0), xp is angle PC'O', and 9 is the 
angle between the line of sight and the motion of the blob. This 
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approximation is always valid because of d > c i 3> r s ^. Here we 
neglect the change of angle O due to the motion of the blob, be- 
cause the distance of the blob is so far away from the observer. 
It is pointed out that we compare the observer time f„b s with the 
time at the center of the blob t' b , i.e., f b s = t' b /T>, where T> is the 
Doppler factor given by T> = 1 /T(l - (3 cos 6), with T being the 
Lorentz factor of the blob motion and /3 = (1 — 1 /r 2 ) 1 / 2 . Then 
we can set up the relation between radiation time t' (time at P) 
and the blob time at the center: 

t> = tL+ d « c '~ d «> *ti + (^.) COS ii>, (7) 

with sufficient accuracy. The retardation is clear in this equa- 
tion. The shell thickness Ar s h increases with the radius of the 
blob, but Ar s h/r s h -C 1 holds for the expansion phase (Sedov 
1969). Thus our treatments are valid. 

The synchrotron emissivity of relativistic electrons is de- 
noted by e v , then the emission from a shell element 
(Ar sh ,A6>',A0') should be AL„ = e ly (r s 2 h Ar sh )sin6' / A6''A0 / , 
where e„(f) = ct,B f^ m " n(j,t)F(x)dx, x = v/uqB^ 2 , F(x) = 
x J x °° K 5 t 3 (z)dz, and K^n is the modified Bessel function 
with order 5/3, vo = 3e/4irm e c = 4.21 x 10 6 , and C3 = 
V3e 3 /(4irm e c 2 ) = 1.87 x 1(T 23 (Pacholczyk 1973). The blob 
is optically thin in our interested bands, which are much higher 
than the synchrotron self-absorption frequency. The emergent 
spectrum from the shell of the expanding blob can be written as 

fix /•27T p'yi 

L(v'X) = c 3 B / smO'de' / d<f>' / H(ON("f,t')F(x)d-y, 

(8) 

where H{() is the Heaviside function: H(Q = for £ < 0, 
whereas H(Q = 1 for £ > 0. Here parameter £ is given by 

C = fb+ — cosV>>0, (9) 
c 

which is incorporated from equations (6) and (7). The intro- 
duction of the Heaviside function is based on the retardation of 
light propagation. Because of the infinitesimal thickness of the 
shell, the retardation of photons to the observer is caused only 
by emission from different points (9', <f>'). It is interesting to 
note that the emission (eq. 8) is independent of the specific ge- 
ometric thickness of the shock. This leads to the avoidance of 
the discussion of the shock thickness. 

We assume the magnetic field changes with radius as 

B{r) = B a {j^j , (10) 

where is the radius of blob at t = t\ . The index n remains 
a free parameter in our model. The flux from the whole blob 
received by a distant observer is obtained by Lorentz transfor- 
mation as 

F(vMs)=^pV'L^,Vt oh ^J , (11) 

where di is the luminosity distance from the source to the ob- 
server. Here we neglect the difference in the directions of emit- 
ting locations in the blob measured from the line of sight. Also 
for simplicity we neglect the difference of Doppler factor at 
different locations in the blob, because the expansion is sub- 
relativistic. 



As a summary of the parameters employed in the present 
model, all of them are listed in Table 1. There are 12 pa- 
rameters describing the processes in the blob. Some of them 
are the initial values of the blob, for example, Z?o, r° h , and /3q. 
With the help of Sedov self-similar solution, these parameters 
at any time will be known. The escape timescale is taken to 
be a free parameter. The dependence of magnetic field on the 
shell radius owing to expansion is an adjustable parameter. The 
timescale of injection ti is a new parameter, and t\ in fact con- 
nects with the initial expansion velocity (3q (eq. 23 below). For 
a non-expanding blob there are 11 parameters (1 - 11). The 
non-expanding blob in fact is homogeneous for a distant ob- 
server, whereas the expanding blob shows the inhomogeneous 
properties as we show below. 

3. QUALITATIVE ANALYSIS 

Generally speaking there are four aspects to describe the 
characteristics of the light curves at different wavebands: 1) 
the symmetry of light curves about the peak flux, 2) the relative 
amplitudes of light curves at the peaks, 3) duration of the bursts, 
which is defined by the e-folding time of vF v , 4) the time lag of 
the light curves, which is defined as the time interval between 
the peak fluxes at different wavebands. Sophisticated observa- 
tions should discover these properties by comparing observed 
spectra with the theoretical models. 

From equation (^) we know that there are four processes 
affecting the observed light curves: blob's expansion which 
increases the blob's size, energy losses of electrons via syn- 
chrotron radiation, photon traveling, and the escape of elec- 
trons from the blob shell. The four processes are competing 
each other, determining the properties of light curves, and are 
hidden in the observed light curves. We first analyze the roles 
of the size, escape, and energy losses in the light curves before 
making detailed numerical calculations. 

3.1. Non-expansion cases 

Considering the complexities of the expanding blobs, we 
first analyze the most simple cases of the non-expansion blobs. 
There are three important timescales, i.e., 1) r syn : energy loss 
timescale due to synchrotron radiation, which is proportional to 
2) ritr ; light-travel timescale through the blob, 3) r esc : the 
timescale of electron escape from the radiation region. Here we 
consider two light curves of photons with frequencies V\ and 
vi (v\ > 1/2). The lifetimes of relativistic electrons which ra- 
diate photons at frequencies is\ and 1/2 are denoted by r s ^' n and 
T S y 2 n , respectively. The total energy radiated at v during lifetime 

r S yn «: 7~' is given by 

AE~ 7^(7,0 cx 7 1_Q . (12) 

For v\ and v%, their total radiated energies are AE\ oc ^f\~ a and 
AE2 fx l\~ a within and r^- n , respectively. Generally the 
size of the blob determines the duration of a flare, if is 
long enough. The escape effect will sharpen the falling pro- 
file and the peak flux, if r esc is short enough. The retardation of 
the peak flux at different bands is obviously determined by the 
above three processes in the blob, i.e., synchrotron loss, light 
travel, and electron escape. The competition among the three 
timescales determines the light curves from the blob. Accord- 
ing to these timescales, we distinguish the following four kinds 
of light curves for a non-expanding blob. 
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(NE1) For the case of r^ n < r^ n < T\ ti < r esc , energy loss 
timescales are shorter than T\ tx and much shorter than r esc . This 
implies that the escape can be neglected. Under such a circum- 
stance, only one process, namely, light traveling, determines 
the observed light curves, because the lifetime of electrons is 
shorter than light travel timescale T\ a . Then the observed light 
curves are mainly due to the light travel effects. The light curves 
should be symmetric about the peak flux, because the rising 
time and the falling time are roughly equal to T\ a . The time 
lag may roughly be zero. The radiation processes of relativistic 
electrons will be hidden. The duration of the flare is roughly 
equal to light travel timescale T\ a -, too. Now we define LCi and 
LC2 as the light curves at frequencies v\ and 1/2, respectively. 
The peak fluxes of LCi and LC2 are denoted by Zq k and ZJf, 
respectively. Thus the ratio of the peak fluxes will be given by 



Affl/Tltr 

L f ~ AE 2 /T ltr 



(l-a)/2 



(13) 



from equation (12). 

(NE2) For the case of r$ n < r ta < r^ n < r esc . The longer 
timescale of escape guarantees the rough conservation of the 
number of radiating electrons in the blob. The light curve at fre- 
quency v\, LCi, is still determined by the blob's geometry, but 
one at v%, LC2, will be determined by radiation processes and 
blob's geometry. The rough symmetry remains in LCi , whereas 
asymmetry in LC2 appears. This is simply owing to the differ- 
ences between rising timescale (due to light travel and energy 
losses) and falling timescale (due to radiation). The duration 
time of LCi will be determined by the light travel timescale 
Tin-, and the duration of LC2 will be modulated by light travel 
timescale and radiation timescale Titc+T^L. The time lag is con- 
trolled by synchrotron radiation and light travel. The duration 
of LCi is still determined by the blob's size, but the duration of 
LC2 will be roughly determined by synchrotron radiation. The 
ratio of the peak flux will be roughly given by 



J 2 



A£i/7ta 



LiT A£ 2 /(Tl tr + T s ? n ) 



1 + - 



Tltr 



(l-a)/2 



(14) 



implying the combined effects of the spectrum, size, and mag- 
netic field. 

(NE3) For the case of T\ a < T S y' n < t%L <C r esc . The shortest 
timescale is the travel timescale, which means that the blob's 
geometry is less important in the light curves. The time lag is 
caused by the differences of energy-loss timescales at v\ and 
V2- The retardation of the flux peak will be of order 

r lag ~ T S y a — T sy ' n . (15) 

The ratio of the amplitudes is determined by the spectrum of 
relativistic electrons, roughly given by 



rPk 

u \ 

rP k 



AE 2 /Tsyn 



(2-a)/2 



(16) 



The duration of a flare is mainly determined by the lifetime of 
electrons radiating at this frequency, which can be roughly ex- 
pressed by t%L/t%* ~ (vi/viY 1 / 1 - The profile of light curves 
may be asymmetric, because the rising timescale is determined 
by the injection, and the falling timescale is determined by en- 
ergy losses. 



(NE4) The case of r s ^ n < r esc < < r hr . The above three 
cases are presumed that the escape timescale is much longer 
than any other timescales. This assumption guarantees the con- 
servation of electron number during the whole lifetime of rela- 
tivistic electrons. The escape effects result in the decrease of the 
electron number (see Eq. ||). Equivalently, the total luminosity 
will be lowered. Now we consider an example of electron es- 
cape. In such a case the escape and light travel effects are im- 
portant in the light curves. If the escape timescale of electrons 
is shorter than light travel and synchrotron radiation timescales, 
the case will be meaningless in observations, because the light 
curves are mainly determined by the escape process. On the 
other hand, if t"1 
plitudes will be given by 



. yn < r esc < 7™, < T\ a holds, the ratio of the am- 



rPk 
rPk 



A£i/7f y ' n 
AE 2 /r esc 



(l-a)/2 



(17) 



which is controlled by the electron spectrum, escape, and mag- 
netic field, and this case will be complicated. The rising profile 
of LC2 is caused by synchrotron radiation, whereas its falling 
profile is by the escape process. The asymmetry appears in light 
curves, and the relative amplitude will be determined by escape 
and radiation processes. 

The above discussions are summarized in Table 2, which 
gives the properties of light curves according to symmetry, du- 
ration, relative amplitude, and time lag. We rule out the case of 
r esc < rmn(ritr, T^,"^), because electrons will directly escape 
without emission. 



3.2. Self-similar expansion 

The expansion of a blob causes the light curves more compli- 
cated because of an additional timescale r exp . The light curves 
from an expanding blob will be given by the above four char- 
acteristic light curves with the modulations by expansion. Gen- 
erally the expansion timescale is longer than that of the light 
travel timescale in the present model, because we are consider- 
ing the non-relativistic expansion. Complexities arise from the 
changes of magnetic field and size of the blob due to expansion. 
Based on the four kinds of light curves described above, it is 
easy to distinguish the new characteristics of light curves from 
an expanding blob. We can compare the expansion timescale 
with the other three timescales; Table 2 is the summary of the 
meaningful light curves according to the four aspects of sym- 
metry, duration, relative amplitude, and time lag. 

The expansion not only enlarges the size of the blob, but also 
weakens the magnetic field. The enlarging size makes the light 
travel timescale longer; the weakening magnetic field dilates 
timescale of synchrotron emission but lowers the output emis- 
sion power. From equation (|l0|), we obtain the synchrotron 
timescale 

2« 



r syn (7) = tL(t) 



(18) 



where T s ° yn = 1.0 x 10 9 /(7#o) s, with r° h and Bq being the initial 
radius and magnetic field of the blob, respectively. For a mono- 
energetic electron with energy 7, the energy loss rate is given 
by 

7(7) = 7o(7) ll) , (19) 
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and the peak frequency of its emission is 

*c(7) = «£(7)(;Sr) ' ( 2 °) 

where 70 = 1.29 x lO" 9 ^ 2 and v c (j) = 4.21 x lOV^o Hz. 
The index n is between 1 and 2 for different models. From this 
equation we see that the synchrotron timescale is very sensitive 
to the expansion law of a blob. This makes the light curves very 
complicated, showing more than one kind of light curves in one 
flare. 

Although the light curves behave similarly to those of the 
non-expanding blob, four timescales such as r syn , Ti tv , r exp , and 
T eS c jointly control the light curves. Comparing them for the 
expanding blob, we obtain seven kinds of light curves and list 
them in Table 2 according to their basic properties such as sym- 
metry, duration, relative amplitude and time lag. 

4. QUANTITATIVE ANALYSIS 

In this section we show our numerical results. As we al- 
ready stressed, the present model includes the expansion effects 
through two processes of enlarging size and weakening mag- 
netic field. We first show the results of a non-expanding blob, 
and then show the results of several cases for different param- 
eters of an expanding blob, and last we compare the results of 
non-expanding and expanding blobs. We list all the cases for 
non-expansion and expansion in Tables 3 and 4, respectively. 
The photon energies are taken to be 0. 1, 1, 10, 100 eV and 1, 2, 
3, 4 keV, which span almost 5 orders of magnitude from near 
infrared, optical, UV, EUV, soft X-ray, to the middle range of 
X-ray bands. The Doppler factor is taken to be V = 10, and 
the total injection energy £0 = 10 45 erg within time interval 
Also 7 m ; n = 10 2 and 7 max = 5 x 10 5 are assumed in all models. 
Finally, Y = 1 /O is assumed. 

4. 1 . Electron energy spectra 

Taking the form of equation (||) as the injection function, 
the electron energy spectrum can be obtained from equation 
(||). Here we should note that the injection function is based 
on Sedov's self-similar solution and the evolution of the elec- 
tron spectrum is self-consistent. Figures 2a and 2b show the 
electron spectra for two cases. For the non-expanding case, the 
number density increases with time for t < t%- The increase is 
linear, because the injection rate is almost constant before time 
?2. Another prominent characteristic is the break at high energy. 
This is due to cooling of synchrotron radiation. The difference 
in the spectral index between the low and high energy is unity. 
This will be reflected in the photon spectral energy spectrum as 
we show in the following figures of photon index evolution. 

There are two possible cases for expansion cases. 1) The in- 
jection is faster than the energy losses. Since the synchrotron 
timescale is dependent on the electron energy, this condition 
will be satisfied for electrons with high energy in the late phase 
of the expansion by considering the magnetic field decays with 
time as t 2 "/ 5 . In such a case the number density will almost lin- 
early increase with time. 2) Synchrotron loss is more rapid than 
expansion (r syn < T exp ). Then the number density of high energy 
electrons will decrease, because the injection is too slow. Such 
a case may appear in the early phase of expansion, because the 
condition (r syn < r exp ) will be satisfied. Then the spectral break 
will appear. In Figure 2b, the above two characteristics are ev- 
ident: in the early phase there is a break due to energy losses, 



whereas in late phase the number density increases linearly and 
the break becomes weak. This will result in different behav- 
iors in the plot of spectrum index versus luminosity in the later 
phase. 

4.2. Non-expanding Blob 

We now consider a non-expanding shocked shell of the blob 
in order to facilitate to compare with the expanding situations, 
although the non-expanding model may not operate in practice. 
Here we assume that the injection still takes the form of equa- 
tion (Q). The size and magnetic field hold constant. Therefore 
nine parameters are in the non-expansion model. Except ti, 
other eight parameters are often used in homogeneous model, 
namely, the magnetic field B, Doppler factor T>, the injected to- 
tal energy Eq, the index of the injected electron spectrum, the 
maximum and minimum Lorentz factors of the electron popu- 
lation, the size of the blob r^ b , and the last is the timescale of 
electron escape from the radiation region. 

In this model the most important timescale is the lifetime of 
electrons radiating synchrotron photons with v. The simple es- 
timation of the lifetime is given by 

r syn = 2.0 x 10 4 B- 3/2 X>- 1/2 ^ 5 1/2 s, (21) 

where all the parameters are in observer's frame, and = 
v/10 15 . The lifetime is insensitive to Doppler factor. We al- 
ways keep the Doppler factor constant in various models. In 
this case the duration time, time lag, and relative amplitude of a 
flare in different bands will be independent of the total injected 
energy. Although there are other seven parameters in the sim- 
ple homogeneous model, the key parameters are r° h , B, and r esc 
which determine the duration time, the time lag, and the rela- 
tive amplitudes in different bands. We show the details in Figs 
3-7. 



4.2.1. Size effects 

Size effects are shown in Figure 3. We change the size as a 
free parameter and fix all other parameters. We take the size of 
the blob, r sh , to be 10 15 , 10 16 ,5 x 10 16 , and 10 17 cm to find the 
effects of the size on light curves. 

Symmetry: The lowest panel in Figure 3 shows the light 
curves for r s h = 10 17 cm. It is interesting to find that the light 
curves in high energy bands (right panel) are symmetric about 
the peak, whereas the rising timescale is much shorter than the 
falling timescale in the low energy band (left panel). This quan- 
titatively confirms the qualitative analysis in §3. The rising and 
falling timescales are determined by the light travel timescale 
ritr as we listed for high energy light curves in Table 2. This 
characteristic can be justified by comparing Tin and r syn . From 
the left panel, we see that the symmetry of light curves is gradu- 
ally broken with the decrease of the size of the blob. We also see 
that the rising timescale decreases rapidly with r S h, because it 
is mainly determined by max(Titr,r syn ). Usually the rising time 
is shorter than the falling timescale for the asymmetric light 
curves. Other asymmetric light curves can be found in Figure 3 
according to the condition listed in Table 2. 

Duration (ra): We find the duration time is very sensitive to 
the size of a blob, because we defined the duration time as the e- 
folding time of vF u . Generally Td will be prolonged by the size. 
The duration time can be roughly given by rj ~ max(T syn ,Titr). 
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The keV photons are radiated from electrons with high energy, 
and it is expected that r syn < T\ a . Thus t$ will be given by T\ a . 
From the figures in the left panel we find that Td of lower energy 
photons is much longer than those in the right panel. The lower 
energy photons produced by electrons with lower energy, and 
it is expected that r syn > T\ tI . Thus the duration time Td will be 
mainly determined by r syn as shown in the left panel, especially 
for the light curves of 0.1, 1, and lOeV photons. Table 2 lists 
the duration properties of LCi and LC2 for the non-expansion 
cases. 

Relative Amplitude: We find the light curves change with r s h 
in several factors. First, for model A4, the condition of NE1 is 
satisfied; all the synchrotron timescales (for high energy light 
curves, see the right panel in Figure 3) are shorter than T\ tr . 
The relative amplitudes are roughly approximated by equation 
(pj|). With the decrease in the blob size, the relative ampli- 
tudes for high energy light curves roughly keep constant but 
they become very complicated for low energy light curves, be- 
cause the condition NE1 is broken. For model A2, r s h = 10 15 cm, 
Titr is shorter than any synchrotron timescales for lower energy 
light curves (left panel in Figure 3), then the condition NE3 is 
satisfied. The relative amplitudes are given by equation (|T6|). 
The upper panel of the left side in Figure 3 indeed shows this 
property. Increasing the size, the relative amplitudes become 
complicated as shown for models A3, Al, and A4. Second, 
the peak flux decreases with the size, if condition NE1 is satis- 
fied. This property can be seen from the right panel of Figure 
3. The left panels show the relative amplitude of peaks at dif- 
ferent bands non-monotonously decreases with the size for low 
energy light curves. This reflects the changes of conditions in 
the blob due to the size of blob. Third, the duration is generally 
prolonged with the increase in the size. All the three features 
can be explained by the fact that the changes of the size of a 
blob modify the light travel time. When the size is large enough 
for a fixed frequency, the broaden light curves will be of small 
amplitude, because the timescale of energy losses of electrons 
is long, which leads to the relative amplitudes of variations at 
lower frequency dominate over high frequency. 

Time Lag: The most prominent properties of these light 
curves are the time lags among different wavebands. From Fig- 
ure 3, we find ri ag is very sensitive to the size with compli- 
cated behaviors. For the symmetric light curves there is zero- 
lag (ri ag = 0), because the rising and falling timescales are given 
by Ti tr as for the case of r s h = 10 17 cm at 1, 2, 3, and 4keV wave- 
bands. However for the asymmetric light curves, the time lag 
becomes complex. As we have already pointed out that the ris- 
ing timescale is roughly determined by max(T syn , ri tr ). Thus for 
high energy light curves the time lag Td is mainly affected by 
the energy losses, if the size is small enough as shown by the 
right panel of Figure 3, which does not show significant time 
lag, because the rising timescale is too short to display by a 
large scale time axis. It is evident to find that time lags increase 
significantly with the size for lower energy light curves from the 
left panel of Figure 3. We still find time lags for r s h = 10 15 cm 
are too small to illustrate in Figure 3. Since r syn > ri tr for 0.1, 
1, and lOeV photons, it is expected that time lag is mainly de- 
termined by the synchrotron energy losses. For intermediate 
cases Td will be jointly given by two processes, i.e., light travel 
and synchrotron radiation. Our numerical results confirm the 
qualitative analysis listed in Table 2. 

Here we have addressed the effects of the blob's size on the 
properties of light curves. These features are key observational 



implications, as we have shown and listed above in detail. Be- 
cause the properties of light curves are highly dependent on 
the observable frequencies, it is desired to compare multiwave- 
length light curves. 



4.2.2. Effects of magnetic fields 

Equations ( 14]), (|l6|), and ( JT^ ) already show the role of mag- 
netic fields. The effects of magnetic field for the non-expansion 
cases are shown in Figure 4. We change the magnetic field as 
a free parameter, while all other parameters are fixed. The role 
of magnetic field can be deduced from the timescale of syn- 
chrotron losses r syn oc Z?" 3//2 [equation (plj)]. The stronger mag- 
netic field, the shorter the timescale of synchrotron losses. This 
leads to 1) the luminosity is enhanced by increasing magnetic 
field, 2) the rising timescale becomes shorter if the size is small 
enough. We fix r s h = 10 16 cm and compare the model Al with 
others to show the effects of magnetic field. Correspondingly, 
the peak flux reaches in advance by increasing the magnetic 
field. This leads to the shortening of time lag T\. dg . This can be 
directly seen in Figure 4. Of course the duration timescale Td 
also shortens by increasing magnetic field. If the magnetic field 
is suitable, symmetric light curves will appear by combined ef- 
fects of synchrotron losses and light travel. For r s h = 10 16 cm, 
we will obtain symmetric light curves, if we increase B until 

T syn < T \ti - 



4.2.3. Effects of electron energy index 

Without a specific acceleration mechanism, we assume that 
the injection spectrum has a power-law with index a. The equa- 
tions in qualitative analysis indicate the importance of index a. 
The roles of the power-law index a of the injected electrons are 
shown in Figure 5. From equation (^]) we know that the lower 
energy electrons are dominant over high energy electrons, if 
a > 2. The power radiated in lower energy will be larger than 
that in high energy bands. This is already shown in Figure 3, 
where we take a = 2.5 in several models. Thus the contribution 
to the rising side of light curves due to the synchrotron radia- 
tion of electrons shifted from higher energy can be neglected, 
because its number is much less than the initial lower energy 
electrons. Time lag is very sensitive to the index a. Lowering 
a leads to the postpone of the flux peak, and shortening time 
lag. The relative amplitude of peak fluxes decreases with the 
increase of a. In such a case the relative amplitude of the peak 
fluxes would be related with the time lag. 



4.2.4. Effects of electron escape 

The role of escape of electrons from the radiation region is 
considered. It is easy to find from equations (|j) and ( |l7| ) that 
the number of electrons decreases with time. With the inclu- 
sion of the size effects, the escape have three factors: the vast 
decreases of fluxes, the shortening of time lag, and the asym- 
metry of light curves. These can be seen in Figure 6, where we 
fix the size of the blob as r s h = 10 I6 cm. The escape effect is 



small on high energy light curves, because t 



syn 



< T„ 



For the 



cases with longer escape timescale, its effects can be neglected. 
For example, T esc /Ti tr = 20 is long enough. Its effects on time 
lag Td for low energy light curves is significant: the decrease of 
Td with increasing T esc as well as the relative amplitude [see eq. 
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4.2.5. Effects of injection timescale 

The effects of injection timescale on the light curves can 
be deduced from equation (Q). The longer the injection, the 
lower amplitudes of peak fluxes as shown in Figure 7, because 
the total injected energy Eq is fixed. The time that fluxes at- 
tain their peak will be postponed with increases of tz- It can 
be found that Td decreases with increasing ti. The symme- 
try strongly depends on the injection timescale ti as shown by 
right-upper panel. The reason is simple: three timescales, t\ v , 
T syn , and ti, are competing in the blob. If ti is short enough 
compared with T syn and T\ tr , the injection can be regarded as im- 
pulsive injection. Then the light curves will be determined by 
the other two timescales r syn and T\ a . For longer duration of 
injection, the light curves will be complicated. For a case of 
?2 > max(T syn ,Ti tr ), the rising timescale will be determined by 
the injection. Generally speaking, the rising timescale will be 
given by max(f2, T syn , T\ a ), whereas the falling timescale will be 
given by max(T syn ,ri tr ) as we already discussed. 



4.2.6. Evolution of photon spectrum 

Figure 8 shows the plot of photon index (a p h) vs. luminosity 
for model Al . The evolution direction is clockwise. We find the 
behaviors of the evolution at 1, 2, 3, and 4keV are quite similar. 
The photon index increases with luminosity (phase 1) before it 
reaches a maximum, and then the index stays constant (phase 
2). Phase 1 is caused by injection and energy loss; the decrease 
of a in phase 1 occurs, because high energy electrons lose en- 
ergy quickly and a break appears in the electron spectrum as 
shown in Figure 2 (left panel), while the increase in the photon 
energy flux is due the injection of electrons. On the other hand, 
in phase 2, the shape of the electron spectrum is steady but the 
total electron number decreases, which leads to the decrease in 
the photon energy flux. It can be found that Aa p h ~ 0.2 at IkeV, 
and Attpj, decreases with increasing photon energy. 

4.3. Expanding Blob 

We have pointed out that the non-expanding model is in fact a 
homogeneous model, because the geometry and magnetic field 
do not change during the evolution of a flare. When a blob ex- 
pands, although the expansion occurs at a constant speed and 
the density is homogeneous everywhere inside the shell in our 
model, the observer will see an inhomogeneous blob, consider- 
ing the photon propagation; the received photons from a blob 
come from different parts of the blob and the different travel 
time from different parts results in emission from an effectively 
inhomogeneous blob. Then the light curves of the expanding 
blob show very complicated behaviors. We refer the expand- 
ing blob as modified homogeneous model. All three timescales 
vitr, Tsyn, and r esc change with time. These result in complicated 
situations for the expanding blob. 

The expansion due to the internal energy can be described 
by Sedov's self-similar solution; the radius of the shocked shell 
reads 

2/5-1 



non-relativistic expansion. This gives 



r sh = r" h max 



1.(1 

h 



(22) 



The parameter t \ is assigned for the validity of the self-similar 
solution. The exact value of this parameter should be obtained 
from solving the dynamics of a blob. Here we simply take it as 
an unimportant parameter in our model. But we constrain it by 



o 

'sh 



5 f3 c 



(23) 



where (3qc is the expansion velocity at t\ . We will carry out the 
research for the expansion with relativistic speed in future. We 
take the index n = 1 , 2 in equation (|To|). The initial speed is set 
by v = ^ h Ai. 

We calculate the models listed in Table 4 in order to find the 
properties of an expanding blob. 



4.3. 1 . Effects of magnetic field index 

The index n of magnetic field will strongly affect the prop- 
erties of light curves. We can understand these effects from 
equations ( ]l8| ) - (p0[). Especially for n = 2, when expansion 
r s h = 2r° h , synchrotron radiation timescale r syn and its mono- 
electron power 7 will be reduced to 1/16 of the non-expansion 
case. It is expected that relative amplitude, time lag, sym- 
metry, and duration time are distorted seriously. Thus the ef- 
fects of expansion will be of significantly observational im- 
plications. Figure 9 shows the light curves of models exp- 
Al and exp-Cl, showing the dependences of light curves on 
the magnetic field index n. In the left panel of Figure 9, it is 
seen that the shapes of light curves are very sensitive to the 
index n. The duration for n = 1 is fundamentally different 
from that for n = 2. Comparing exp-Cl and exp-Al with non- 
expanding Al models (the parameters are the same besides ex- 
pansion), we can find interesting differences. Between Al and 
exp-Cl, we find four properties as follows. 1) the luminos- 
ity is significantly lowered by expansion, because the magnetic 
field weakens via expansion as B oc r~". 2) the relative ampli- 
tudes are changed, L pk (100eV)/L pk (10eV) > 1 for expansion, 
whereas L pk (100eV)/L pk (10eV) < 1 for non-expansion case. 
The similar properties can be found by comparing L pk (100eV) 
and L pk (leV) etc. 3) Duration time Td is changed significantly. 
Model exp-Cl has much longer Td than Al, because the syn- 
chrotron energy loss rate become smaller for weakening mag- 
netic field. 4) Time lag Ti ag is shortened via expansion, because 
the flux peak time occurs in advance. This effect is compli- 
cated, because the rising timescale will be controlled by three 
timescales such as r syn , T\ a , and r exp . For high energy light 



curves, r syn will be shorter than T\ tI and T exp . Thus the ris- 
ing timescale will be mainly determined by r syn as compar- 
ing with non-expansion cases. Also time lag Ti ag is expected 
[it is not significant for keV light curves, but we can compare 
L pk (0.1keV) and L pk (lkeV)]. It should be noted that h has a 
strong effect on the peak flux time. But it is beyond the scope 
of this paper to discuss the value of t \ , because we do not solve 
hydrodynamics of the expansion. 

Between Al and exp-Al, light curves are almost changed 
completely by the expansion with n = 2. First, output luminos- 
ity weakens. 2) The relative amplitudes are changed similar to 
exp-Cl. 3) Duration time Td is completely different between 
Al and exp-Cl. The luminosity weakening is due to the the 
rapidly weakening magnetic field rather than the energy losses 
of electrons. Thus the duration will be controlled by expan- 
sion. The falling profiles are t~'"(m « 0.6 ~ 0.7) rather than 
linear decrease, namely they become dim with powerlaw of 
time. From equation dlq), we know the synchrotron timescale 
becomes very long (T syn oc r^ h ). This leads to the decrease in 
luminosity owing to the decrease of magnetic field rather than 
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the energy loss of electrons. Thus this is completely different 
from the other cases for n = 1. Under such a circumstance, 
we quantitatively estimate the light curves from equation (|T^); 

vF v oc AE/r exp , then 

vF v oc jy(i-«)/2 f (i-")«/5 ; (24) 

where we use ^ oc Bj 2 oc f ~ ln l i ^ 1 . We see the falling profile t~° 6 
(for n = 2 and a = 2.5) and the frequency-dependence i/ 1-Q V 2 
This crude approximation is in agreement with the detail calcu- 
lations shown in Figures 9 and 10 for the model exp-Al, exp- 
Cl, and exp-C2. 4) There is almost no significant time lag be- 
tween the light curves both for high and low energy light curves. 
This is because r exp <C r syn . 5) The peaks are very sharp. This 
is a new characteristics caused by index n. Expansion leads to 
strong n-dependence of synchrotron timescale. After t = t\, the 
decline of light curves is caused by weakening magnetic field 
rather than energy losses of relativistic electrons. 

4.3.2. Effects of initial expansion speed 

We fix the initial radius of a blob r° h , but adjust t\ in order to 
test the effects of initial speed. Figure 1 1 shows the light curves 
for different velocity /3 = 0.22c (t\ = 6 X 10 5 s) and /3 = 0.17c 
(h = 8 X 10 5 s). 

Generally the expansion enlarges the size of a blob, and thus 
leads to similarity with the cases of larger radius if the expan- 
sion is fast enough. Figure 10 clearly shows the effects of initial 
expansion speed. The effects are clearer for high energy light 
curves (keV bands). The lower initial velocity (3q, the faster 
the peak flux is attained. The peak flux significantly lowers, 
because the faster expansion leads to rapid decrease in mag- 
netic field. Then the duration time ta for faster initial expansion 
speed is shorter than that for slower initial expansion speed. On 
the other hand, time lag is not changed much by initial expan- 
sion speed. Time lag ri ag is jointly controlled by r syn , T] tr , and 
T exp . We summarize in Table 2. For high energy light curves 
the rising timescale may be mainly determined by r syn , and the 
falling timescale is controlled by T\ tI and T exp . Thus the sym- 
metry of light curves will be generally broken. The expand- 
ing blob will always have asymmetric light curves, because the 
timescales are changing during the expansion. 

4.4. Evolution of spectrum 

Figure 11 illustrates the trajectory in the photon index vs. 
flux plane for an expanding blob. The upper panel shows the 
evolution of the spectrum for n = 2. It can be seen that the 
behaviors are quite different from those of the non-expansion 
cases. The spectral evolution is similar for non-expanding and 
expanding blobs, before the luminosity attains its maximum. 
The photon index a p h increases with increasing luminosity. 
However, the evolution behaviors are different after then. We 
see aph continues to increase with decrease in luminosity for an 
expansion blob. The lower panel shows the comparison for the 
cases of n = 1 and n = 2. We find that a p h becomes larger values 
for n = 2. We also find that the behavior difference becomes 
large with decreases of photon energy. These results are mainly 
because of the decrease in magnetic field. 

5. CONCLUSIONS 

We presented a very simple Sedov-expansion model of a 
coasting blob and studied the received spectrum radiated from 



it. The multiwaveband light curves were calculated for the 
lower energies such as 0.1, 1, 10, and 100 eV, and higher en- 
ergies such as 1,2, 3, and 4 keV. Main four processes deter- 
mine the observed light curves, namely, the evolution of the 
energy spectrum of relativistic electrons, the size (initial and 
enlarging size) of the radiation region, bulk motion of the blob, 
and electron escape from the radiation region. According to the 
timescales, r syn , Ti tr , r exp , and r exp , we find 1 1 kinds of light 
curves, which are listed in Table 2 with a concise summary of 
their properties including symmetry, duration, relative ampli- 
tude, and time lag. To facilitate to compare our model with 
observations, we stress the followings: 

(1) In X-ray band the rising profile of light curves is deter- 
mined by the timescale of light travel in a blob, and the declin- 
ing profile is symmetric, if the timescale of synchrotron radia- 
tion is short. 

(2) In EUV band the rising timescale of the light curve is the 
sum of Ti tr , r sy n V , and the lifetime of electrons radiating EUV 
photons. The declining profile is determined by r sy ^ v . 

(3) The observed UV light curve profile is determined by the 
whole processes, and the rising timescale is roughly the sum of 
Tin-, T ' S yn V ' anc l r syn- The declining timescale is determined by 
the minimum of r esc and rHJ'. 

(4) Light curves from an expanding blob depend on how 
magnetic field decreases with expansion. When the decrease 
of luminosity is steep, it may be because of the rapid decrease 
of magnetic field and it is possible to estimate the index n of 
magnetic field. 

(5) The relationship between the duration time, rising time, 
decline time, and the relative amplitudes in different wavebands 
is expected to be found from the multiwaveband observations. 
It is useful to test the undergoing flare mechanisms by using the 
model presented here. 

(6) For the rapid decays of magnetic field, the falling pro- 
files can be described as vF v oc ^( 1_a )/ 2 f( 1_Q; )"/ 5 ? which is sig- 
nificantly different from the slowly decaying one without ex- 
pansion. Thus the profile can be used to test the magnetic field 
in the blob. 

These properties apply for the emission spectrum which has 
the synchrotron tail in X-ray bands; also note that this situation 
depends on the value of 7 max and V, too. 

In this paper we calculated the light curves of an expanding 
blob interacting with its surroundings. The electron injection 
function was obtained analytically and self-consistently based 
on the Sedov's self-similar expansion. We neglected the re- 
acceleration of electrons which lose energy via synchrotron ra- 
diation. Since the energy loss timescale may be comparable 
with the expansion timescale (or the injection timescale), the 
acceleration region and the radiation region in fact are almost in 
the same region. Then the re-acceleration of electrons is possi- 
ble (Schlickeiser 1984). The pile-up energy will shift with time. 
We will study this subject in future. The time-dependent injec- 
tion with a temporal profile may lead to more complex light 
curves. It is needed to study the detailed acceleration process 
to interpret the time-dependent injection pattern, which is far 
beyond the scope of the present paper. In our present calcula- 
tions we used the same approximation made by Dermer (1998) 
and Georganopoulos and Marscher (1998) who neglected in- 
verse Compton scattering. The validity of this approximation 
has been discussed in Georganopoulos and Marscher (1998). If 
the expansion is relativistic then the difference of Doppler fac- 
tor at different locations in a blob can not be neglected. We will 
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Table 1 . The input parameters in the present model 



1. Bq blob's initial magnetic field 

2. V blob's Doppler factor 

3. Eq Energy of relativistic electrons injected into a blob 

4. a index of relativistic electrons at injection time 

5. 7 max maximum energy of electrons in a blob 

6. 7 m i n minimum energy of electrons in a blob 

7. r° h initial radius of shocked shell of a blob 

8. T esc timescale of electron escape from the radiation region 

9. t\ timescale of expansion before self-similar solution is applied 

10. ?2 timescale of injection 

11. f3o initial expansion velocity of a blob in units of c 

12. n index of magnetic filed 



Table 2 Summary of kinds of light curves at two frequencies v\ and vi (nu < t„ p ) 









symmetry 






duration 


















non-exp 




exp 


non-exp 


exp 


relative amplitude 


time 


ag 


non-expansion 


expansion 


LCI 


LC2 


LCI 


LC2 


LCI 


LC2 


LCI 


LC2 


non-exp 


exp 


non-exp 


exp 


(1) r,« < T«, < Tte < TescCNEl) 


Tsyn < Ts'yn < Utr < T CX p < T CSC (E1) 


yes 

(Tltr,Tltr) 


yes 

(TUr,Tltr) 


yes 

(T|tr,T| tr ) 


yes 

(TltP,Tto) 


ltr 


ltr 


ltr 


ltr 


sp+ltr 


sp+ltr 


no 


no 


(2) < IK < r«, < r csc (NE2) ■ 


T„y„ < Th, < 7%, < Tcxp < T«c(E2) 
<y'„ < Tin < T cxp < T«, < T CSC (E3) 


yes 

(Tltt.Tltr) 


no 

(Tltt.T,^,) 


yes 

(Tltr,Tl tr ) 

yes 

(T|tr,Tl tr ) 


no 

(Tl tt ,T s «) 

no 

(Tcxp,^^) 


ltr 


ltr + syn 


ltr 
ltr 


syn 
exp 


sp + ltr 


ltr + syn 
exp + syn 


syn 


syn + ltr 
syn + exp 


(3) r ltr < r s " y '„ < r^, < r csc (NE3) ■ 


' T,,, < T$, < T,« < T cxp < T CSC (E4) 
Tto < T^! n < Tcxp < T^, < T es c(E5) 
. T|„ < Tcxp < 7$, < T^, < Tcsc(E6) 


no 


no 

(t" j t" 2 ) 

v 'syn! ' syn' 


no 

(Tltr.T^) 

no 

(Tltr.T^) 

no 

(T e xp,T s ^,) 


no 

^ 'syni ' syn' 

no 

(T e xp,T s ^ n ) 

no 

(T" : T n 1 
V 'syn) ' syn' 


t" 1 

' syn 


t" 1 + t" 2 

' syn 1 ' syn 


syn 
syn 
exp 


syn 
exp 
exp 


sp 


sp 

syn + exp 
exp + syn 


syn 


syn 
syn + exp 
syn + exp 


(4)r s ^<r eS c<[7l tt ,r s «](NE4) 


T s " y „ < T„c < [Tltt.T^J < Tcxp(E7) 


no 

(TbscT^) 


no 

(Tltr,T es c) 


no 

(Tfy'n.Tcc) 


no 

(T!tr,T e sc) 


ltr+syn 


ltr+esc 


ltr+esc 


ltr+esc 


sp+ltr+esc 


esc+syn 


syn+ltr+esc 


esc+syn 



Note: "ltr" indicates the light travel effect, "exp" does the blob's expansion, "syn" does the radiation of electron due to synchrotron, "sp" does the spectrum, and "esc" does the escape from the blob. The symbol "+" means that the oberved quanty is jointly 
controlled by the several processes inidcated by the abbreviated words above. In the "symmetry" column we also indicate the rising timescale (left) and falling timescale (right) below the "yes" or "no". 



Table 3. Parameter values for non-expansion models 



Model Al A2 A3 A4 Bl B2 B3 B4 CI C2 Dl D2 El E2 



a 


2.5 


2.5 


1.5 2 


2.5 


f 2 (xl0 6 ) s 


5 


5 




2.5 1 


B G 


0.1 


0.2 0.3 0.4 0.6 0.1 






r° h (xl0 16 )cm 


1 0.1 5 


10 1 1 






7~esc/ Tltr 


10 


5 20 


10 





Note: The default values of the parameters are referred to the same as their left values in this table. 



Table 4. Parameters for expansion models 



Model 


exp-Al 


exp-A2 exp-Cl exp-C2 


a 


2.5 




f 2 (xl0 6 ) s 


5 




Bo G 


0.1 




n 


2 


1 


r° h (xl0 16 ) cm 


1 


1 




10 




fi(xl0 5 ) s 


6 


3 6 8 



Note: The default values of the parameters are referred to the same as their left values in this table. 
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Z 

A 




Fig. 1 . — The illustration of the expanding blob in its coasting frame. The black shell is the shock formed via the collision between 
the expanding blob and the inter-clouds media in broad line region. The shocked shell is very thin compared with the radius of the 
blob. The photons arrive at the observer at different time due to the different distances. Here 9 is the angle between the line of sight 
and the blob motion. 
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Fig. 2. — The evolution of the electron energy spectrum. Left panel (Fig. 2a) is for a case of non-expansion (model Al), whereas 
the right panel (Fig. 2b) is for an expansion case (model exp-Al). For model Al, a broken power law appears after a certain time 
interval which is determined by timescales of energy loss and escape. For model exp-Al, the break energy is higher than that of 
model Al, because the energy loss is slower owing to the decreasing magnetic field. 
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Fig. 3. — The light curves for non-expansion blob for different values of r s h = 10 , 10 16 , 5 x 10 , and 10 17 cm from top to 
bottom. The values of parameters for the models are listed in Table 3. 
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The light curves for different magnetic field for a non-expanding blob: model Al, Bl, and B2 from top to bottom. 
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Fig. 5. — The light curves for different values of the index of the injected electron spectrum for a non-expanding blob: model Dl, 
D2, and Al from top to bottom. 
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The light curves for different escape timescales for a non-expanding blob: model CI, Al, and C2 from top to bottom. 
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Fig. 7. — The light curves for different injection timescales characterized by t2 for a non-expanding blob: model E2, El, and Al 
from top to bottom. 
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Fig. 8. — The trajectory in the photon index a and luminosity plane for a non-expanding model Al. The evolution is clockwise. 
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Fig. 9. — The light curves for expansion cases: exp-Cl (top) exp-Al (bottom). In the top panel, a non-expanding model Al is also 
shown by thin curves for comparison. 
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Fig. 10. — The light curves for expansion cases with exp-Cl (t\ = 6 x 10 5 s) and exp-C2 {t\ = 8 x 10 5 s). Thick curves are for exp-Cl 
and thin curves are for exp-C2 
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